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ABSTRACT 

We describe calculations of the evolution of an ensemble of small 
planetesimals in the outer solar system. In a solar nebula with a mass of several 
times the Minimum Mass Solar Nebula, objects with radii of 100-1000 km 
can form on timescales of 10-100 Myr. Model luminosity functions derived 
from these calculations agree with current observations of bodies beyond the 
orbit of Neptune (Kuiper Belt objects). New surveys with current and planned 
instruments can place better constraints on the mass and dynamics of the solar 
nebula by measuring the luminosity function at red magnitudes, itir > 28. 



Subject headings: solar system: formation - Kuiper Belt 



1. INTRODUCTION 



Several remarkable discoveries have renewed interest in solar system formation. Recent 
surveys have detected many small icy bodies beyond the orbit of Neptune (e.g., [Jewitt & 



Luu 1993 ; Williams et al. 1995| ; Jewitt et al. 199^ ; Luu et al. 1997 ; [Gladman fc KavelaarsI 



|1997| ; [Chiang fc Brown 1999|) . Assuming a geometric albedo of 4%, these Kuiper Belt 
objects (KBOs) have radii of 50-400 km; the derived size distribution implies a significant 
population of smaller objects. With semi-major axes of 40-50 AU and orbital inclinations 
of 0°-30°, the orbits of known KBOs suggest an annulus of planetesimals formed in situ 
and left over from the planetary formation epoch (|Holman fc Wisdom 1993]) . The presumed 
structure of this annulus resembles the dusty disks recently discovered around several 
nearby stars ( [Smith fc 'lerrile 1984| ; [Aumann et al. 1984| ; [J ayawardhana et at. 1998| ; [Koerner 
\et al. 1998^ [Greaves et al. 1998| ). Planets similar to those in our solar system have not 
been detected in any of these disks, but direct images and radial velocity measurements of 
other nearby stars already imply the existence of more than one dozen extra-solar planets 
of several Jupiter masses ( [Latham et al. 1989[ ; [Marcy fc Butler 1996[ ; [Cochran et al. 1997[ ; 
Noyes et al. 1997[ ; [DelFosse et al. 1998[ ; for a review see [Marcy 1999[ ) . 



These discoveries challenge planet formation theories. Most theories presume that 
planets grow by accretion of small planetesimals in a gaseous circumstellar disk ( [Safronov 
[1969[ ; [Goldreich fc Ward 1973| ; see also [Lissauer fc Stewart 1993[ , [Boss 1997[ , and references 
therein). Hitherto, numerical studies have focused on the formation of the prototypical 
terrestrial and gas giant planets. Earth and Jupiter ([Greenberg et al. 1978| , 1984; [Nakagawa 
et al. 19831 ; phtsuki et al. 1988[ ; [Wetherill fc Stewart 1989| , 1993; [Barge fc Pellat 1990[ 



Ruden fc Pollack 199 1[ ; [Pollack et al. 199q ; [Weidenschilling et al. 1997[) . If the initial disk 
mass is comparable to the Minimum Mass Solar NebulaQ, these calculations often have 
difficulty producing objects similar to the known terrestrial or gas giant planets during 
the estimated disk lifetime of ~ 10-30 Myr ( Pr'oUack et al. 1996[ ; [Weidenschilling et al. 



[1997[ ). This problem is exacerbated in the outer solar system, where numerical calculations 
yield formation times exceeding 100 Myr for 500-1000 km radius KBOs ( Stern 1995 , 1996; 
[Stern fc Colwell 1997a] ,b). KBOs must form on shorter timescales in parallel with Neptune. 
Otherwise, Neptune's gravity increases the velocities of nearby planetesimals, including 
those in the inner Kuiper Belt, on timescales of ~ 10 Myr (jMalhotra 1996[) . This process 



^ The Minimum Mass Solar Nebula has a surface density E = S]o(^/^o) '^^^, where Eo is the surface 
density of solid material at i?o = 1 AU. We adopt Sq = 45 g cm^^ ( Hayashi 1981 ; see also Weidenschilling 
1977; Bailey 1994). This definition yields a total mass of solids Mq « 100 Me inside the orbit of Neptune 
{R < 30 AU) and Mq « 10 Me in the inner part of the Kuiper Belt {R = 32-38 AU), where 1 M^ = 6 x 10^^ 
g is the mass of the Earth. 



prevents the growth of KBOs with radii exceeding 100-200 km, because large velocities 
hinder agglomeration. 

We recently began to consider KBO formation in the outer solar system using an 
evolution code that follows planetesimal growth in the annulus of a circumstellar disk. 
Initial results indicate that KBOs can form at 30-50 AU on timescales of 10-100 Myr 
in disks with 1-3 times the Minimum Mass Solar Nebula when collisional disruption of 
planetesimals is unimportant (Kenyon & Luu 1998; hereafter |KL98|) . Further calculations 
with an algorithm that includes disruptive processes lead to similarly short timescales for 
a wide range of initial conditions (Kenyon & Luu 1999; hereafter P:<.L99| ). Here, we briefly 
summarize these new results, compare the theoretical model with current observations, and 
make predictions for comparison with future observations of KBOs. 



2. MODEL 

Our accretion code is based on the particle-in-a-box method, where planetesimals are a 
statistical ensemble of bodies with a distribution of horizontal and vertical velocities about 
Keplerian orbits ( ^afronov 1969 ). We perform calculations for a single annulus of width Aa 



centered at a heliocentric distance a. We approximate the continuous distribution of particle 
masses with i discrete batches having particle populations riiit) and total masses Mi{t). 
The horizontal and vertical velocity dispersions are hi{t) and Vi{t) ( [Wetherill fc Stewart] 



1993| ). The average mass of a batch, mi{t) = Mi{t)/ni(t), changes with time as collisions 
add and remove bodies from the batch. This procedure conserves mass and provides a 
statistical method to follow the growth of ^ 10^° small planetesimals into a few planets. 
Detailed n-body calculations confirm the basic features of particle-in-a-box calculations for 
the early stages of planet growth described here ( |lda fc Makino 19^ ; [Kokubo fc Ida 19"^ ). 



To evolve the initial size distribution in time, we calculate collision rates for the 
coagulation equation, determine the outcome of each collision, and compute velocity 
changes due to collisions and long-range gravitational interactions (see KL99). Each 
two-body collision can produce (1) merger into a single body with no escaping debris (very 
low impact velocity), (2) merger into a single body with escaping debris ('cratering'; low 
impact velocity), (3) rebound with or without cratering (modest impact velocity), or (4) 
catastrophic disruption into numerous smaller bodies (high impact velocity). The collision 
outcomes depend on the ratio of the impact energy Qf to the disruption energy Qd of two 
colliding planetesimals ( [Greenberg et al. 1978| ; [Wetherill fc Stewart 1993| ; Pavis et al. 1994] ) . 



Collisions with Qf > Q^ disrupt planetesimals into many small fragments. Collisions with 
Qf < Qd yield a merged planetesimal and some small fragments if the collision velocity V^ 



exceeds the miniinuin velocity for cratering Vj. Collisions with Qf < Q^ and Vc < Vf yield 
a merged planetesimal with no cratering debris. We use an energy- scaling formalism to 
compute Qd as the sum of the intrinsic tensile strength 5*0 and the gravitational binding 
energy ( pavis et at. 198^ , 1994). The intrinsic strength is the dominant component of Qd 
for bodies with ri ^1 km; gravitational binding dominates 5*0 for larger bodies. 

For each collision, a velocity evolution algorithm distributes the kinetic energy among 
the resulting bodies and then accounts for collisional damping, kinetic energy transfer 
during elastic collisions ("dynamical friction"), angular momentum transfer during elastic 
collisions ("viscous stirring"), and gas drag ( |Hornung et al. 1985| ; see also |Wetherill fc 



Stewart 1993| ; |KL98| ). Dynamical friction tries to enforce equipartition of kinetic energy 



between mass batches; viscous stirring increases the velocities of all bodies. Gas drag 
removes objects from the annulus and reduces the velocities of small objects which are well 
coupled to gas in the disk. 



We tested the code against analytical solutions and published numerical results (|KL98 



P<.L99| ). We reproduced previous calculations for accretion at 1 AU ( |Wetherill fc Stewart 



1993| ) and collisional disruption of pre-existing large KBOs at 40 AU ( pavis fc 1^'arinelTa 



1997 ). Our calculations match analytical solutions well when the mass spacing between 
successive batches, S = mi^i/rrii = 1.1-1.4. Numerical solutions lag the analytic results by 
~ 10% when 6 = 1.4-2. The timescale to produce objects of a given size increases with 6, 



because poorer resolution prevents growth of large objects (see [KL98|) . 



Table 1 lists basic input parameters. The input cumulative size distribution Nc has the 
form Nc oc rf°, with initial radii r^ = 1-80 m. The total mass in the annulus is Mq; Mq ~ 
10 Me for a Minimum Mass Solar Nebula. All batches start with the same initial velocity. 
We tested a range of initial velocities corresponding to initial eccentricities of cq = 10~^ 
to 10~^, as is expected for planetesimals in the early solar nebula ( [Malhotra 1995| ). The 
adopted mass density, po = 1-5 g cm~^, is appropriate for icy bodies with a small rocky 
component. The fragmentation parameters - Vf, Sq, Qc, fxE, Ci, and C2 - are adopted from 
earlier work. KL99 describe these parameters in more detail. 

To provide observational constraints on the models, we note that the known Kuiper 
Belt population contains at least one body with a radius of ~ 1000 km (Pluto), and ~ 10^ 
KBOs with radii Tj ^ 50 km between 30-50 AU. The cumulative size distribution of known 



KBOs can be fitted with A^^ oc r'-"'"', with qobs = 3 ± 0.5 ( |Jewitt et al. 1998t see also phiang 
|fc Brown 1999| ). Successful models should reproduce these observations on timescales 
comparable to (a) the estimated lifetimes of the solar nebula and gaseous disks surrounding 



nearby young stars, <^ 10^ yr (p^ussell et al. 1996| ; Partmann et al. 199^ ) and (b) the 



formation timescale for Neptune, <^ 10* yr ( [Lissauer et al. 1996|) . 



3. NUMERICAL RESULTS 

We separate the growth of KBOs into three regimes. Early in the evolution, frequent 
collisions damp the velocity dispersions of small bodies. These bodies slowly grow into 1 km 
objects on a timescale that is approximated by ri km ~ 8 Myr (Mq/IOMe) (eo/lO"^)"'^^. 
This linear growth phase ends when the gravitational range of the largest objects exceeds 
their geometric cross-section. This "gravitational focusing" enhances the collision rate by 
factors of 10-1000. The largest objects then begin a period of "runaway growth", when 
their radii grow from ~ 1 km to ;^ 100 km in several Myr. During this phase, dynamical 
friction and viscous stirring increase the velocity dispersions of the smallest bodies from ~ 
1 m s~^ up to ~ 40 m s~^. This velocity evolution reduces gravitational focusing factors 
and ends runaway growth. The largest objects then grow slowly to 1000+ km sizes on 
timescales that again depend on Mq and Cq. Column (5) in Table 2 lists timescales to form 
Pluto-size objects Tp as a function of the input parameters Mq, 6, Cq, and go- 
Fig. 1 shows cumulative size distributions for a model with Mq = 10 Me, qo = 3, 
and 5*0 = 2 X 10^ erg g~^. The shapes of these curves depend on two competing physical 
processes: (1) growth by mergers and (2) erosion by high velocity collisions. In this example, 
collisions result in growth because the velocity dispersion is less than the catastrophic 
disruption threshold. However, the collision velocity exceeds the cratering threshold Vf. 
Cratering adds debris to all low mass batches. Gas drag removes material from low mass 
batches (rj ^ 10 m), but is ineffective at removing larger objects. The size distribution thus 
becomes shallower at small masses. At large masses, mergers produce a group of growing 
planetesimals with a steep size distribution. Once gravitational focusing becomes effective, 
the largest of these objects 'run away' from the rest of the ensemble to produce a smooth 
power law with a maximum radius r^ax- As the evolution proceeds, r^ax increases but the 
slope of the smooth power law remains nearly constant. 

The main features of these results depend little on the input parameters. All 
calculations produce two cumulative power law size distributions connected by a transition 
region having an 'excess' of planetesimals (the "bump" in the curves in Fig. 1). The 
characteristic radius of this transition region increases from 0.3 km at Cq = 10""^ to 3 km at 
Co = 10^^. If fitted with a power law of the form Nc oc r^ '^^ at small masses, the cumulative 
size distribution follows the predicted limit for collisional evolution, qj = 2.5 ( pohnanyi 



1969| ). We perform least-square fits to obtain qf at larger masses; column (6) of Table 2 
lists derived values for qf along with the la error. Column (7) lists the radius range for 
each fit. The results are surprisingly independent of the input parameters. We find the 
small range qf = 2.75-3.25 for calculations with Mq = 1-100 Me, Cq = 10~^ to 10~^, go = 
1.5-4.5, and 5*0 = 10 erg g~^ to 3 x 10^ erg g~^. This model result is consistent with the 



observed slope, qobs = 3 ± 0.5 (e.g., Jewitt et al. 1998). 



4. COMPARISONS WITH OBSERVATIONS 

As shown in Table 2, several calculations meet the success criteria defined in §2. 
Annuli with Mq ^ IOMe produce Pluto-sized objects on short timescales, Tp ^ 50 Myr (for 
6o ^ 10^'^). Models with smaller initial masses or larger initial eccentricities form Plutos on 
longer timescales, Tp ^ 50 Myr. Plausible ranges of other input parameters - such as qo, 
Sq, and fxE ^ yield ±20% variations about these timescales. The results are insensitive to 
Vf and other collision parameters ( [KL99[ ). 



The crosses in Figure 1 compare our calculations directly with several observational 
constraints. The cross at Tj = 50 km indicates the number of KBOs with Tj ^ 50 km 
estimated from recent ground-based surveys (|Jewitt et al. 199^ ; see also [Chiang fc Brown 



1999| ); the one at rj = 10 km shows limits derived from a single, controversial measurement 



with Hubble Space Telescope {HST; [Cochran et al. 1995| , 1998; [Brown et al. 1997D . The 



third cross plots limits at r^ = 1 km derived from theoretical attempts to explain the 
frequency of short-period comets from the Kuiper Belt ( Davis fc Farinella 19971 Duncan 



[fc Levison 1997| ; [Levison fc Duncan 1997[ ). Our predictions agree with ground-based 



surveys at 50 km and theoretical limits at 1 km, but fall a factor of ~ 10 short of the HST 
measurement at 10 km. 

To compare with observations in more detail, we predict the luminosity function 
(LF) of KBOs directly from the computed number distribution. We use a Monte Carlo 
calculation of objects selected randomly from the cumulative size distribution Nc. We 
assign each object a distance from the Sun (Iq and a random phase angle (3 between the 
line-of-sight from the Earth to the object and the line-of-sight from the Sun to the object. 
This phase angle lies between 0° and a maximum phase angle that is distance-dependent. 
The distance of the object from the Earth is then cIe = dQCOsP — (1 + dQ{cos'^P — 1))^^^. 
We derive the red magnitude of this object from a two parameter magnitude relation for 
asteroids, mp^KBO = Ro + 2.5 log (^1/^2) — 5 log tkbo, where Rq is the zero point of the 
magnitude scale, tkbo is the radius of the KBO, ti = 2dQdE, and ^2 = i^((l — g)4>i + 5^02) 
( [Bowell et al. 195P| ). In this last expression, u is the albedo, and g is the slope parameter; 
01 and 02 are phase functions that describe the visibility of the illuminated hemisphere of 
the object as a function of /5. We adopt standard values, u = 0.04 and g = 0.15, appropriate 
for comet nuclei ([Jewitt et al. 199^ ). The zero point Rq is the apparent red magnitude of 
the Sun, rripQ = —27.11, with a correction for the V-R color of a KBO, Rq = rn^Q + 
(5(V-R)ft:BO- Observations suggest that KBOs have colors that range from roughly —0.1 
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to 0.3 mag redder than the Sun. We treat this uncertainty by allowing the color to vary 
randomly in this range. 

The important parameters in the model LF are the distributions of input sizes (derived 
from the accretion calculations), distances, and orbital parameters. We assume KBOs are 
evenly distributed between "Plutinos," objects in 3:2 orbital resonance with Neptune having 
semimajor axes of 39.4±0.2 AU, and "classical" KBOs with semimajor axes between 42-50 
AU. The distance parameters are set by observations ( |Jewitt et al. 199^ ). This distance 



distribution is different from the 32-38 AU adopted for the coagulation calculations. Several 
tests show that accretion results at 42-50 AU are identical to those at 32-38 AU, except 
that the timescale to produce Pluto-sized objects is 50%-100% longer. To compute the 
model LF from the Monte Carlo magnitude distribution of KBOs, we scale the mass in 
the 32-38 AU annulus to match the mass in a 42-50 AU annulus, add in an equal number 
of Plutinos, and divide by the sky area. The distribution of KBO orbital parameters is 
poorly known. We adopt circular orbits to derive magnitudes; the LF is insensitive to other 
choices. We assume orbital inclinations of i = 0° to 5° to compute the sky area, which is a 
compromise between the i ~ 0°-5° of classical KBOs and the i ~ 10°-30° of Plutinos. The 
model LFs scale inversely with sin i. 

Figure 2 compares several models with the observed LF. The left panel shows models 
with Co = 10~^ and different masses; the right panel shows models with the mass of a 
Minimum Mass Solar Nebula and different Cq. Model LFs with the Minimum Mass and 
any cq agree with current observations. The good agreement of all models at tur < 20, 
where the uncertainties are largest, depends on the assumed maximum radius in the model 
distribution. We picked 1000 km for convenience. Model LFs for rriR > 20 are independent 
of this choice. 

To quantify the comparison between models and observations, we fit model LFs to 
log Tj{mR) = ainiR — mo) over 20.5 < niR < 26.5. Table 3 lists the fitted a and mo as 
a function of the mass in classical KBOs (in units of the Minimum Mass Solar Nebula), 
Co, the inner annulus boundary Rin, and the outer annulus boundary Rout- The small 
range in a for model LFs agrees with published values derived from observation^. The 
model a is independent of the relative numbers of Plutinos and classical KBOs, and the 
distance distribution of classical KBOs. The observed zero-point of the LF, mo ~ 23.2-23.5, 



^Gladman et al. (1998) report a = 0.761q '^5 and ttiq = 23.4^Q'^g from a maximum likelihood analysis of 
previous surveys with magnitude limits, 20 < mn < 28. Surface densities in their Table 3 yield a « 0.6 and 
mo « 22.4. Jewitt et al. (1998; see also Luu & Jewitt 1998) quote a = 0.54 ± 0.04 and mo = 23.2 ± 0.10 for 
20 < m^ < 26. Chiang & Brown (1999) prefer a = 0.52 ± 0.02 and mo == 23.5 for 20 < ttir < 27; they note 
that the slope depends on which survey data are used in the fit. 



favors models with masses comparable to the Minimum Mass Solar Nebula and any initial 
eccentricity. These data rule out models with < 30% of the Minimum Mass at the 3a level. 
Smaller Plutino fractions require larger masses: if Plutinos are 10%-25% of the total KBO 
population, as indicated by recent observations ( Jewitt et al. 190^ ), the needed mass is 2-4 
times the Minimum Mass. 

There are two main uncertainties in comparing our model LFs with the data, the 
evolution of the KBO LF with time and the current orbital parameters of KBOs. The 
initial mass in KBOs was larger than implied by a direct comparison between the data and 
model LFs, because large velocity collisions and dynamical encounters with Neptune have 
eroded the Kuiper Belt over time ( [Holman fc Wisdom 1993| ; pavis fc 1^'arinella 1997) ; see 



also [Levison fc Duncan 1993| ; puncan et al. 1995| ). Erosion from high velocity collisions 
probably does not change the slope of the LF significantly. Massive KBOs with r^ ^ 50 
km {rrifi <^ 26-27) are probably safe from coUisional disruption ( Davis fc Farinella 1997 ). 



Disruption of smaller bodies depends on the unknown bulk properties and the poorly 
known orbital parameters of KBOs. These uncertainties are not important for comparisons 
of models and observations for niji ^ 26-27, but can bias future comparisons at fainter 
magnitude limits. Gravitational perturbations from Neptune should affect all KBO masses 
equally and simply reduce the total mass in KBOs with time ( [Holman fc Wisdom 1993| ). 
Despite the uncertainty in the total amount of mass lost from the Kuiper Belt, we are 
encouraged that the mass needed to explain current observations of KBOs is at least the 
Minimum Mass Solar Nebula. Future calculations will allow us to place better constraints 
on the initial mass in the Kuiper Belt. 

The uncertain distribution of KBO orbital parameters also affects the initial mass 
estimates. Our assumption of KBOs uniformly distributed in distance dQ, orbital 
eccentricity e, and inclination i is probably incorrect for Plutinos in specific orbital 
resonances with Neptune. Larger adopted volumes for current Plutinos require larger initial 
disk masses in the Kuiper Belt. A uniform distribution is probably reasonable for classical 
KBOs, but the observed range in (Iq and i is not well-known. Allowing classical KBOs to 
occupy a larger range in semi-major axis reduces our mass estimates; a larger range in sin 
i increases our mass estimates. We suspect that the uncertainties currently are a factor of 
~ 2-3. Future large area surveys will provide better knowledge of KBO orbital parameters 
and allow more accurate models for the observed LF. 

In addition to the reasonably good fit for 20.5 < rriji < 26.5, our calculations predict 
1-5 'Plutos' with rriR < 20 over the entire sky. This number is uncertain, because we do not 
understand completely the mechanism that ends accretion and sets the maximum size of 
KBOs. Our calculations indicate that planetary accretion at 35-50 AU is self-limiting: once 



objects reach radii of ~ 1000 km, they stir up smaller bodies sufficiently to limit additional 
growth. The formation of nearby Neptune should have also limited the growth of the largest 
bodies ( |lVlorbidelli fc Valsecchi llJ^ ). Better constraints on the radial distribution of 500+ 
km KBOs would test the relative importance of these two mechanisms. 

Observations at fainter magnitude limits will provide additional constraints on KBO 
formation. Imaging data acquired at the Keck and Palomar telescopes detect KBOs with 
rriij ~ 25-26.5, where models with cq ^ 10~^ predict the LF to rise sharply. The apparent 
lack of a significant upturn in the LF at mji < 25 implies cq ^ a few x 10~^. In contrast, 
the current limit on the KBO population at rriji > 28 implies a substantial population of 10 
km radius KBOs which is inconsistent with our calculations. Deeper ground-based surveys 
or new HST data could resolve the controversy surrounding this observation and place 
better constraints on cq. Finally, the proposed Next Generation Space Telescope {NGST) 
will probe the size distribution of 1 km radius KBOs where models with cq ~ 10^'^ predict 
a sharp upturn in the observed LF. If such small bodies can survive for the age of the solar 
system, NGST observations would provide important constraints on the initial mass and 
dynamics of the outer solar system. 



We thank B. Bromley for making it possible to run our code on the JPL Cray T3D 
'Cosmos' and the HP Exemplar 'Neptune' and for a generous allotment of computer time 
through funding from the NASA Offices of Mission to Planet Earth, Aeronautics, and Space 
Science. Comments from F. Franklin, M. Geller, and M. Holman greatly improved our 
presentation. 



REFERENCES 

Aumann, H. H. et al. , 1984, ApJ, 278, L23 

Bailey, M. 1994, In Asteroids, Comets, Meteors 1993, edited by A. Milani, M. DiMartino, 
and A. Cellino, Kluwer Academic Publishers, Dordrecht, p. 443 

Barge, P., and Pellat, R. 1990, Icarus, 85, 481 

Boss, A. P. 1997, Science, 276, 1836 

Bowell, E., Hapke, B., Domingue, D., Lunime, K., Peltoniemi, J., and Harris, A. W. 1989, 
in Asteroids II, edited by R. P. Binzel, T. Gehrels, and M. S. Matthews, Tucson, 
Univ. of Arizona Press, p. 524. 

Brown, M. E., Kulkarni, S. R., & Liggett, T. J. 1997, ApJ, 490, L119 

Chiang, E. I., & Brown, M. E. 1999, AJ, 118, No. 3 (|astro-ph/9905292D 



Cochran, A., Levison, H., Stern, S. A., & Duncan, M. 1995, ApJ, 455, 352 

Cochran, A. L., Levison, H. F., Tamblyn, P., Stern, S. A., & Duncan, M. 1998, ApJ, 455, 
L89 

Cochran, W. D., Hatzes, A. P., Butler, P. R., & Marcy, G. W. 1997, ApJ, 483, 457 

Davis, D. R., Chapman, C. R., Weidenschilling, S. J., & Greenberg, R. 1985, Icarus, 62, 30 

Davis, D. R., & Farinella, P. 1997. Icarus, 125, 50 

Davis, D. R., Ryan, E. V., & Farinella, P. 1994, Planet. Space Sci., 42, 599 

DelFosse, X., et al. 1998, A&A, 338, L67 

Dohnanyi, J. W. 1969, J. Geophys. Res., 74, 2531 

Duncan, M. J., & Levison, H. F. 1997, Science, 276, 1670 

Duncan, M. J., Levison, H. F., & Budd, S. M. 1995, AJ, 110, 3073 

Gladman, B., & Kavelaars, J. J. 1997, A&A, 317, L35 

Gladman, B., Kavelaars, J. J., Nicholson, P. D., Loredo, T. J., & Burns, J. A. 1998, AJ, 
116, 2042 

Goldreich, P., & Ward, W. R. 1973, ApJ, 183, 1051 

Greaves, J. S. et al. 1998, ApJ, 506, L133 

Greenberg, R., Wacker, J. F., Hartmann, W. K., & Chapman, C. R. 1978, Icarus, 35, 1 

10 



Greenberg, R., Weidenschilling, S. J., Chapman, C. R., Davis, D. R. 1984, Icarus, 59, 87 

Hartmann, L., Calvet, N., Gullbring, E., & D'Alessio, R 1998, ApJ, 495, 385 

Hayashi, C. 1981, Prog Theor Phys Suppl, 70, 35 

Holman, M. J., & Wisdom, J. 1993, AJ, 105, 1987 

Hornung, P., Pellat, R., & Barge, P. 1985 Icarus, 64, 2951 

Ida, S., & Makino, J. 1992, Icarus, 96, 107 

Irwin, M., Tremaine, S., & Zytkow, A. N. 1995, AJ, 10, 3082 

Jayawardhana, R. et al. 1998, ApJ, 503, L79 

Jewitt, D., & Luu, J. 1993, Nature, 362, 730 

Jewitt, D., Luu, J., & Chen, J. 1996, AJ, 112, 1225 

Jewitt, D., Luu, J. X., & Trujillo, C. 1998, AJ, 115, 2125 

Kenyon, S. J., & Luu, J. X. 1998, AJ, 115, 2136 (KL98) 

Kenyon, S. J., & Luu, J. X. 1999, AJ, 118, No. 2 ( |astro-ph/9904TT5| ) 



Koerner, D. W., Ressler, M. E., Werner, M. W., & Backman, D. E. 1998, ApJ, 503, L83 

Kokubo, E., & Ida, S. 1996, Icarus, 131, 171 

Kowal, C. T. 1989, Icarus, 77, 118 

Latham, D. W., Stefanik, R. P., Mazeh, T., Mayor, M., & Murki, G. 1989, ApJ, 339, 38 

Levison, H. P., & Duncan, M. J. 1990, AJ, 100, 1669 

Levison, H. P., & Duncan, M. J. 1993, ApJ, 406, L35 

Levison, H. P., & Duncan, M. J. 1997, Icarus, 127, 13 

Lissauer, J. J., Pollack, J. B., Wetherill, G. W., & Stevenson, D. J. 1996, in Neptune and 
Triton, edited by D. P. Cruikshank, M. S. Matthews, & A. M. Schumann, U. of 
Arizona Press, Tucson, 37 

Lissauer, J. J., & Stewart, G. R. 1993, In Protostars and Planets III, edited by E. H. Levy 
and J. I. Lunine, U. of Arizona Press, Tucson, p. 1061 

Luu, J. X., & Jewitt, D. 1988, AJ, 85, 1256 

Luu, J. X., & Jewitt, D. 1998, ApJ, 502, L91 

11 



Luu, J. X., Marsden, B., Jewitt, D., Trujillo, C. A., Hergenother, C. W., Chen, J., & Offutt, 
W. B. 1997, Nature, 387, 573 

Malhotra, R. 1995, A J, 110, 420 

Malhotra, R. 1996, AJ, 111, 504 

Marcy, G. W. 1999, in The Physics of Star Formation and Early Stellar Evolution, edited 
by C J. Lada and N. Kylafis, Dordrecht, Kluwer, in press 

Marcy, G. W., & Butler, R R. 1996, ApJ, 464, L133 

Morbidelh, A., & Valsecchi, G. B. 1997, Icarus, 128, 464 

Nakagawa, Y., Hayashi, C., Nakazawa, K. 1983, Icarus, 54, 361 

Noyes, R. W. et al. 1997, ApJ, 483, Llll 

Ohtsuki, K., Nakagawa, Y., Nakazawa, K. 1988, Icarus, 75, 552 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., & Greenzweig, 
Y. 1996, Icarus, 124, 62 

Ruden, S. P., & Pollack, J. B. 1991, ApJ, 375, 740 

Russell, S. S., Srinivasan, G., Huss, G. R., Wasserburg, G. J., & Macpherson, G. J. 1996, 
Science, 273, 757 

Safronov, V. S. 1969, Evolution of the Protoplanetary Cloud and Formation of the Earth 
and Planets, Nauka, Moscow [Translation 1972, NASA TT F-677] 

Smith, B. A., & Terrile, R. J. 1984, Science, 226, 1421 

Stern, S. A. 1995, AJ, 110, 856 

Stern, S. A. 1996, AJ, 112, 1203 

Stern, S. A., & Colwell, J. E. 1997a, AJ, 114, 841 

Stern, S. A., & Colwell, J. E. 1997b, ApJ, 490, 879 

Weidenschilling, S. J. 1977, Ap Sp Sci, 51, 153 

Weidenschilling, S. J., Spaute, D., Davis, D. R., Marzari, P., Ohtsuki, K. 1997, Icarus, 128, 
429 

Wetherill, G. W., & Stewart, G. R. 1989, Icarus, 77, 300 

Wetherill, G. W., & Stewart, G. R. 1993, Icarus, 106, 190 



12 



Williams, I. P., O'Ceallaigli, D. P., Fitzsimmons, A., Marsden, B. G. 1995, Icarus, 116, 180 



This preprint was prepared with the AAS IfTjTiX macros v4.0. 

13 



Table 1. Basic Model Parameters 



Parameter 


Symbol 


Value 


Width of annulus 


6a 


6 AU 


Initial velocity 


Vo 


0.45-45 m s"^ 


Particle mass density 


Po 


1.5 g cm~3 


Relative gas velocity 


V 


30 m s"^ 


Time step 


St 


1-250 yr 


Number of mass bins 


N 


64-256 


Mass spacing of bins 


5 


1.25-2.0 


Minimum velocity for cratering 


Vf 


1 cm s~^ 


Impact strength 


So 


2 X 10^ erg g"^ 


Crushing energy 


Qc 


5 X 10'' erg g~^ 


Fraction of KE in ejecta 


Ike 


0.05 


Coefficient of restitution 


Cl 


10-2 


Coefficient of restitution 


C2 


10-3 



Table 2. Model Results at 32-38 AU 



Mo{Me) 


6 


eo 


qo 


Tp (Myr) 


Qf 


Range (km) 


1 


1.4 


10-^ 


-3.0 


448 


-2.86 ±0.04 


2-300 


10 


1.4 


10-^ 


-3.0 


20 


-2.90 ±0.02 


1-930 


1 


1.4 


10-3 


-3.0 


893 


-2.76 ±0.05 


6-400 


3 


1.4 


10-3 


-3.0 


184 


-2.78 ±0.03 


5-600 


10 


1.4 


10-3 


-3.0 


37 


-2.91 ±0.03 


6-800 


30 


1.4 


10-3 


-3.0 


10 


-3.02 ±0.03 


6-930 


100 


1.4 


10-3 


-3.0 


3 


-2.97 ±0.03 


7-600 


10 


1.4 


10-2 


-3.0 


428 


-3.15 ±0.10 


50-700 


100 


1.4 


10-2 


-3.0 


25 


-3.23 ±0.06 


20-600 


10 


1.25 


10-3 


-1.5 


40 


-2.97 ±0.02 


7-700 


10 


1.25 


10-3 


-3.0 


35 


-3.03 ±0.03 


9-650 


10 


1.25 


10-3 


-4.5 


30 


-2.90 ±0.02 


4-750 



Table 3. Luminosity Function Parameters 



Mo/M. 



MMSN 



eo 



R,n (AU) Rout (AU) 



a 



mo 



0.3 
1.0 
1.0 
1.0 
3.0 



10-3 

10-^ 
10-3 

10-2 

10-3 



42 
42 
42 
42 
42 



50 
50 
50 
50 
50 



0.56 ± 0.01 24.03 ± 0.16 

0.58 ± 0.01 23.36 ± 0.12 

0.57 ± 0.02 23.16 ± 0.15 

0.60 ± 0.01 23.53 ± 0.18 

0.58 ± 0.01 22.42 ± 0.17 



0.3 
1.0 
1.0 
1.0 
3.0 



10-3 
10-4 
10-3 
10-2 
10-3 



42 
42 
42 
42 
42 



60 
60 
60 
60 
60 



0.56 ± 0.01 23.50 ± 0.13 

0.57 ± 0.01 22.85 ± 0.11 

0.56 ± 0.02 22.63 ± 0.13 

0.63 ± 0.02 23.31 ± 0.19 

0.58 ± 0.01 21.92 ± 0.17 
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Fig. 1. — Cumulative size distributions as a function of time for a model with Mq = IOMe 
and Co = lO""^. The evolution time for each curve is listed in the legend. Crosses indicate 
observational and theoretical constraints on the size distribution at radii of 50 km, 10 km, 
and 1 km as described in the text. 
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Fig. 2. — Comparison of model luminosity functions of KBOs with observations. Data are as 
indicated in the legend of each panel. The open circle with the central dot is the position of 
Pluto for an adopted albedo of 4%; other observations are from Cochran et al. (1998; HST), 
Irwin et al. (1995; 195), Kowal 1989 (1989; K89), Luu & Jewitt (1988; LJ88), Gladman et 
al. (1998; G98), Luu & Jewitt (1998; LJ98), Jewitt et al. (1998; JLT98) and Chiang & 
Brown (1999; CB99). Error bars for each datum - typically a factor of 2-3 - and the upper 
limit from Levison & Duncan (1990) are not shown for clarity. The lines plot luminosity 
functions for models with (a) left panel: cq = 10^^ and Mq ^ 0.3 (dot-dashed), 1.0 (solid), 
and 3.0 (dashed) times the Minimum Mass Solar Nebula and (b) right panel: a Minimum 
Mass Solar Nebula with cq = 10~^ (dashed), eo = 10"^ (sohd), and cq = 10~^ (dot-dashed). 
A Minimum Mass Solar Nebula has Mq ~ 12 Me within R = 42-50 AU. The pair of vertical 
solid hues indicates the planned magnitude range accessible to NGST. 



